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Abstract: We construct rotating black hole solutions in Einstein-Gauss-Bonnet theory in 
five spacetime dimensions. These black holes are asymptotically flat, and possess a regular 
horizon of spherical topology and two equal-magnitude angular momenta associated with 
two distinct planes of rotation. The action and global charges of the solutions are obtained 
by using the quasilocal formalism with boundary counterterms generalized for the case of 
Einstein-Gauss-Bonnet theory. We discuss the general properties of these black holes and 
study their dependence on the Gauss-Bonnet coupling constant a. We argue that most of 
the properties of the configurations are not affected by the higher derivative terms. For 
fixed a the set of black hole solutions terminates at an extremal black hole with a regular 
horizon, where the Hawking temperature vanishes and the angular momenta attain their 
extremal values. The domain of existence of regular black hole solutions is studied. The 
near horizon geometry of the extremal solutions is determined by employing the entropy 
function formalism. 
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1. Introduction 

In d > 5 dimensions, the gravity action may be modified to include higher order curvature 
terms while keeping the equations of motion to second order, provided the higher order 
terms appear in specific combinations In five dimensions, this leads to the so-called 
Einstein-Gauss-Bonnet (EGB) theory, which contains quadratic powers of the curvature. 
The Gauss-Bonnet (GB) term appears in the one-loop corrected effective action of heterotic 
string theory [0, B, and also in the low-energy effective action for the compactification of 
M-theory on a Calabi-Yau threefold Q. Inclusion of this term in the action leads to a 
variety of new features (see the recent reviews |5|, ^]). In particular, the black holes of EGB 
theory do not in general obey the Bekenstein-Hawking area law, but the entropy formula 
includes a new contribution coming from the higher curvature terms in the action j?], ||. 
Also, for d = 5, the GB term implies the existence of a branch of small static black holes 
which are thermodynamically stable ||. 

However, the complexity of EGB theory makes the task of finding closed form solutions 
a highly non-trivial task. The only known analytic solutions correspond to the counterparts 
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of the Schwarzschild-Tangherlini black holes [10, 11], and a variety of physically interesting 



solutions were found only numerically, (e.g. the black strings [12, [L3|] and the black rings 
fl4|| ). In particular, no rotating EGB solutions are known yet in closed form, and it was 
proven in [15] that the Kerr-Schild ansatz may not work in Lovelock theory. Nevertheless, 
a number of partial results support the idea that EGB generalizations of the Myers- Perry 
(MP) solutions [16] actually exist [17, 18]. 



The main purpose of this work is to present numerical evidence for the existence of a 
class of asymptotically flat, rotating solutions in d = 4+1 EGB theory. These solutions have 
a spherical horizon topology and two equal-magnitude angular momenta. This restriction 
leads to a system of coupled nonlinear ordinary differential equations, which are solved 
numerically within a nonperturbative approach. The same approach has been employed 
to construct Einstein-Maxwell ]i~9|, |20[], Einstein-Maxwell-Chern-Simons [21], and Einstein- 



Yang-Mills [22 1 rotating black hole solutions in higher dimensions. Also, the existence of 



the AdS counterparts of the solutions addressed in this work was established in [23]. 

Our results suggest that the EGB rotating black holes with two equal angular momenta 
share most of the features of the corresponding Einstein gravity MP solutions. In particular, 
for a fixed value of the GB coupling constant, the sets of black hole solutions terminate 
when they reach an extremal black hole, where the Hawking temperature vanishes and the 
equal-magnitude angular momenta assume their extremal values. These extremal black 
holes obey the attractor mechanism assuming their near horizon geometry involves an 
AdS2 factor. As a new feature, a GB term in the action leads to the existence of a branch 
of rotating black holes with a positive specific heat at constant angular velocity at the 
horizon. 

This paper is organized as follows: in the next Section we recall the EGB action 
and present the ansatz and boundary conditions. Section 3 contains a discussion of the 
counterterm approach for asymptotically flat solutions in EGB theory, together with the 
computation of the action and global charges for rotating black holes with equal magnitude 
angular momenta. We present our numerical results in Section 3, and exhibit the physical 
properties of these solutions and their domain of existence. The near horizon geometry of 
the extremal solutions is also studied there by using the entropy function formalism 
We conclude in Section 4 with some further remarks. 

2. The model 
2.1 The action 

We consider the Einstein-Hilbert action supplemented by the GB term: 



where R is the Ricci scalar and 



L GB = R 2 - m^R^ + R^ T R^ T (2.2) 



- 2 - 



is the GB term. Variation of the action with respect to the metric tensor results in 
the equations of EGB gravity 



R*» - 2 R 9^ + \ R w = > ( 2 - 3 ) 



where 



H,_ IV — 2(R fJ-t7KT R I/ e7KT — 2R fipucr R pcr — 2R fJ/7 R a l/ + RR^ V ) - -LcbQ^u • (2.4) 

For a well-defined variational principle, one has to supplement the action (^[l]) with the 
Gibbons-Hawking surface term 

m _ L_ f * n 

IdM 



H ' = "g^ / d*xV-hK , (2.5) 



and its counterpart for EGB gravity || 

where h ab is the induced metric on the boundary, K is the trace of the extrinsic curvature 
of the boundary, G ab is the Einstein tensor of the metric h a b and J is the trace of the 
tensor 

J ab = - {2KK ac K c h + K cd K cd K ab - 2K ac K cd K db - K 2 K ab ) . (2.7) 

An interesting feature of EGB gravity is the presence of two branches of static solutions, 
distinguished by their behaviour for a — > [pi]. In this paper we shall restrict our analysis 
of rotating EGB black hole solutions to those, whose static limit corresponds to the branch 
of static solutions with a well defined Einstein gavity limit. 

2.2 The ansatz and boundary conditions 

While rotating EGB black holes will generically possess two independent angular momenta 
and a more general topology of the event horizon 1 we restrict here to configurations with 
equal-magnitude angular momenta and a spherical horizon topology. A suitable metric 
ansatz reads [19j 

dr 2 

ds 2 = — — + g(r)d6 2 + Mr) sin 2 9 (dtpi - w(r)dt) 2 + h(r) cos 2 6 (d<p 2 - w(r)dt) 2 (2.8) 
fir) 

+ (g(r) - h(r)) sin 2 6 cos 2 9{d^ - dip 2 ) 2 - b(r)dt 2 , 

where 9 E [0, 7r/2], (991,992) £ [0, 2tt], and r and t denote the radial and time coordinate, 
respectively. Also, this metric ansatz admits a simpler expression in terms of the left- 
invariant 1-forms <Tj on S 3 , with 



dr 2 1 1 1 

ds 2 = — + -g(r)(a 2 + a 2 ) + -h(r){a 3 + 2w(r)dt) 2 - b(r)dt 2 , (2.9) 



1 Recently, also static black rings were obtained in d — 5 EGB gravity by employing a numerical approach 
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where o\ = cos ipd9 + sin ij) sin 9d(f>, 02 = — sin ipd9 + cos ip sin 9d<fi, u 3 = cfa/> + cos 9d<f> and 

20 = 9^ x -4>i = <P, 0i + <p2 = 1>. 

For such solutions the isometry group is enhanced from R x U(l) 2 to R x U(2), where R 
denotes the time translation. This symmetry enhancement allows to factorize the angular 
dependence and thus leads to ordinary differential equations. 

Without fixing a metric gauge, a straightforward computation leads to the following 
reduced action for the system 



i-cS 



a 



drdt L eff , with L cS = Le + -^ l gb, 



(2.10) 



with 



b g V 9 



Lgb = \I J -^-[ —b'g' + 2(4g - 3/i)(^f- + hw' 2 ) - ^-b'h'g' 2 



T 



'2„../2 



2h 



-^fhg w 



(2.11) 



where a prime denotes a derivative with respect to r. The corresponding equations for the 
metric functions b, f, h and w are found by taking the variation of A e g with respect to 
a, b, f and g and by fixing afterwards the metric gauge. (This procedure is equivalent to 
solving the EGB equations directly, but it is technically simpler.) The equation for w{r) 
admits the first integral 



/, ffhf a fg' 2 



12/t., 
16 H )) = const., 

g 



(2.12) 



which is useful in the numerical calculations. The remaining equations are rather long, and 
we do not include them here. In the following 2 we fix the metric gauge by taking g(r) = r 2 . 

Also, the EGB equations (|2.3| ) imply the following relations which will be important 
in the later discussion 



1 



sin 



1 



—^(4 + ^/4) 

: J_ [LA 

2r 2 V bh dr 



(2.13) 




hw' 



r 2 + a(f-4+ 3 L 



rv 1 
R\ + -(Hi + -L GB ) 




r 2 (hww' — b') 



A- LA 

2r 2 V bh dr t 

+ «((/- 4 +^ + ^r) ft +(4-72 



(2.14) 



f)hww' + rfhw 



The horizon of these black hole solutions is a squashed S 3 sphere. It resides at the 
constant value of the radial coordinate r = r# > 0, and is characterized by f(rji) = 



2 The numerical calculations are performed in two gauges. Besides the gauge g(r) = r 2 , we also employ 
the gauge choice in jH], ^| (with b(r) = f(r), f = f(r)/m(r), g(r) = r 2 m/f(r), h(r) — nr 2 /f(r), and 
w(r) = w(r)/r), corresponding to isotropic coordinates. 
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fr(rff) = 0. Restricting to nonextremal solutions, the following expansion holds near the 
event horizon: 

f(r)= fi(r-r H )+0(r-r H ) 2 , h(r) = h H + 0(r - r H ), (2.15) 
b(r) = b\(r — r#) + 0(r — rjj-) 2 , w(r) = VLh + w\(r — ru) + 0(r — rjj) 2 . 

For a given event horizon radius, the essential parameters characterizing the event horizon 
are /i, b\, hH, &h and w\ (with fi > 0, b\ > 0), which fix all higher order coefficients 
in ( 2.15 ). (These constants are related in a complicated way to the global charges of the 
solutions.) The (constant) horizon angular velocity Qh is defined in terms of the Killing 
vector x = d/dt + Q\d/d(p\ + 0,^8 fdipi which is null at the horizon. For the solutions 
within the ansatz ( |2.8| ), the horizon angular velocities are equal, f2i = Q% = ^H- 

A straightforward calculation gives the following asymptotic expansion for the metric 
functions, involving three parameters U, V and W 

6(r) = 1 + y 2 + ( 2 ^ 2 " WV + 3^ 2 | ) ±g + (U 2 V +U - 2W 2 a) ^ + 0(±), 



+ (§^(W 2 +UV) + ^(W 2 - WV)^ 1 + O(^), 



h(r) = r 2 + ^- + (^(W 2 +^V) + ^(W 2 - aHVja) £ (2.16) 

+ ^ f (5V - 12W 2 )(W 2 + UV) - 7U(W 2 - 2UV)aj ^ + O(^) , 

w w(v-^f) Ln( 1 x 

w ( r ) = 7? ^ + (^To)' 

which guarantees that the Minkowski spacetime background is approached at infinity. 

2.3 Known solutions and slowly rotating black holes 

For the metric ansatz (|2.8|), the EGB field equations (|2.3|) possess two well known exact 



solutions. First, the MP black holes [16] with equal-magnitude angular momenta are found 
for a = (i.e., without GB term). Expressed in terms of the event horizon radius and the 
horizon angular velocity 3 (which are the control parameters in our numerical approach), 
this solution reads 

f(r) = 1 - \ 2 + rl f H 2 ( rjL X , Kr) = r 2 ( 1 + (^Y ^ 2 

J V ' 1 - r 2 tt 2 \ r ) 1 - r 2 n 2 \ r J W I V r J 1 - r 2 fi£ 



Therefore, for a MP black hole, the relevant parameters in the event horizon expansion 
( |2~T5|) are 

h = 2(1 " 'I oil ^ = A (l-2« »i = ^(rX-l), (2-18) 
r H (l-r H n H ) r H r H 



3 This metric is usually expressed in terms of the mass parameters M and the angular momentum 
parameter a, with M = r|j/2(l — r 2 H Q, 2 H ) and a = r 2 H Q.u- 



- 5 - 



while the constants U, V and W in the far field expansion ( |2,16| ) have the following ex- 
pression 4 

r 6 O 2 r 2 r 4 Ou- 

V = 1 r _ H , H o2 , ^ = ~ 1 _ r f o2 , W = 1 y 2 ^ 2 . (2.19) 

The d = 5 MP black holes with equal-magnitude angular momenta emerge smoothly 
from the static Schwarzschild-Tangherlini black hole when the event horizon velocity is 
increased from zero. For a given event horizon radius, the solutions exist up to a maximal 
value of the horizon angular velocity = l/y/2rj{. Expressed in terms of the mass- 
energy E and the equal-magnitude angular momenta \ J\\ = IJ2I = J, this bound reads 
27irJ 2 /8G < E 2 . The extremal solution saturating this bound has a regular but degenerate 
horizon. Further details on the properties of MP black holes with equal-magnitude angular 
momenta are found in 25 , pBf . 



The second exact solution corresponds to the generalization |1C] of the static Schwarzschild 
solution with a GB term and reads (note our restriction to the branch with Einstein gravity 
limit) 

w (r) = 0, h(r) = r 2 and f(r) = b(r) = 1 + T — (l - I + ^{a + 2r 2 H )^ . (2.20) 

This solution exists for all ru > and a > —rjj (with ru the event horizon radius). 
Without entering into details, let us mention the existence of some substantial differences 
between the thermodynamics of the static d = 5 EGB black hole solutions 5 and their 
Einstein gravity counterparts (restricting to the physical case a > 0). If the black holes 
are large enough, r 2 u > a, then they behave like their Einstein gravity counterparts since 
the specific heat is negative, C p = Tn{dS / 8Th) < 0. A different behaviour is found for 
rjj/a < 1, since in that case Th — th and thus C p > 0. This implies the existence of a 
branch of small five-dimensional EGB black holes that is thermodynamically stable Q (see 
e.g. H for a review of these aspects). 

To the best of our knowledge, there is no exact solution of the EGB equations de- 
scribing non-static configurations 6 . However, slowly rotating black holes can be found by 
considering perturbation theory around the static solution ( |2.20| ) in terms of the rotation 



parameter a (see e.g. [17|). For our case the slowly rotating solution then contains in its 



non-diagonal metric elements the function w(r) 



a(2r 2 TT + a) 

w(r) = \ H ' =^ - ; (2.21) 

a(2r%+a) 



r 4 1 + \ 1 + 



that is linear in the perturbative parameter a, while the other metric functions remain 
unchanged to this order in a. 



4 Note that since V = — W 2 /U, there are only two free parameters in the far field expansion, which fix 
the mass and equal-magnitude angular momenta of the solutions. 

For d > 5, the thermodynamical properties of the solutions are similar to the Einstein gravity case. 

6 See however, the results in JTHJ . However, the solution there is very special (e.g. the value of the GB 
coupling constant is fixed by the cosmological constant) and does not describe a black object. 
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3. The boundary counterterm approach in d = 5 Einstein-Gauss-Bonnet 
theory and the global charges 



3.1 A counterterm for the d = 5 asymptotically flat space 

It is well known that the gravitational action contains divergences even at the tree-level - 
that arise from integrating over the infinite volume of spacetime. A common approach - 
background subtraction - uses a second, reference spacetime to identify divergences which 
should be subtracted from the action. After subtracting the (divergent) action of the 
reference background, the resulting action will be finite. 

At a conceptual level, the background subtraction method is not entirely satisfactory, 
since it relies on the introduction of a spacetime which is auxiliary to the problem. In some 
cases the choice of reference spacetime is ambiguous - for example NUT-charged solutions 
(see e.g. the discussion in |27j). For asymptotically AdS spacetimes, this problem is solved 
by adding additional surface terms to the action |28|] . These counterterms are built up with 
curvature invariants of a boundary dM (which is sent to infinity after the integration) , and 
thus they do not alter the bulk equations of motion. This yields a finite action and mass 
of the system. 

The generalization of this procedure to the asymptotically flat case was considered in 
1 29 , 30, 31]. Moreover, as discussed in p2[, a renormalized stress-tensor can be defined by 



varying the total action (including the counterterms) with respect to the boundary metric. 
The conserved quantities can be constructed from this stress tensor via the algorithm of 



Brown and York [33 



However, all studies in the literature of asymptotically flat configurations have con- 
sidered the case of counterterms in Einstein gravity only 7 . We have found that for five- 
dimensional asymptotically flat EGB solutions with a boundary topology S 3 x R, the action 
can be regularized by the following counterterm: 

1 



let 



(fi\/-W(K), (3.1) 

8M 



8ttG 

where 7Z is the Ricci scalar of the induced metric on the boundary hij and 



y(K) = J-K(l + -K). (3.2) 



2 V 9 

For a = 0, this reduces to the counterterm expression proposed in |31] for Einstein gravity 
with the same asymptotic structure. 

Varying the total action (which contains the Gibbons-Hawking boundary term) with 
respect to the boundary metric hij, we obtain the divergence-free boundary stress-tensor 

T ab = T^i^r = s^G (K ab - h ab K + f (Qrf - \Qh ab ) - 2^TZ ab + 2^h ab - 2h ab U^> + * ;afe ) , (3.3) 

where K ab is the extrinsic curvature of the boundary and [Q] 

Qab = 2KK ac K c b - 2K ac K cd K db + K ab (K cd K cd — K 2 ) + 2KK ab + KK ab - 2K cd lZ cadb - ATZ ac K c b , 

7 The general expression of the counterterms and the boundary stress tensor in EGB theory with a 
cosmological constant is presented in [p5|, |36|. 
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while 7Z a bcd, Ttab denote the Riemann and Ricci tensors of the boundary metric. 

Provided the boundary geometry has an isometry generated by a Killing vector a 
conserved charge 



£2 ? = j> d 6 S l (3.4) 



can be associated with a closed surface £ [27]. Physically this means that a collection of 



observers on the hypersurface with metric hij all observe the same value of 0e, provided 
this surface has an isometry generated by £. 



To test the counterterm (3.1), we have verified that the known expressions for the mass 
and action are recovered for the case of the Schwarzschild-Tangherlini solution in EGB. 
Neglecting the second term in ( |3.2| ) (i.e. considering only the Einstein gravity countert- 
erm) results in the occurrence of an unphysical constant term proportional with a in the 
expressions of I c i and E. 

However, one should remark that the proposal ( |3.1| ) inherits all ambiguities already 
present in the Einstein gravity case. In particular, there is no rigorous justification for 
considering that expression, and, in fact the results in |3]J show that the counterterm choice 
in Einstein gravity is not unique. This suggests that other more complicated expressions 
are possible for the EGB case as well 8 . Moreover, unlike the AdS case, the counterterm 
depends on the boundary topology (for example, we have verified that the coefficient in 
front of ([3.1|) is different for black string solutions). 



These problems are avoided in the proposal put forward by Mann and Marolf in [37] for 



asymptotically flat spacetimes and further generalized in [38] to extended objects, like black 
strings and p— brane spacetimes. In this approach, the conserved quantities are constructed 
essentially using the electric part of the Weyl conformal tensor. The generalization of the 
Mann and Marolf prescription to Lovelock gravity solutions is an interesting open problem. 

3.2 The action and global charges of solutions 

The computation of the boundary stress-tensor T a f, for the solutions of interest in this 
work is straightforward, and, for the asymptotic expression (2.16), we find the following 
expressions for the relevant components 

^^^W>> ^^W), zt-±£ + oww> 

The energy E and the equal-magnitude angular momenta J are the charges associated with 
the Killing vectors d/dt, d/d<pi, and d/d(f2, respectively. Computed according to ( |3~3| ) 
these quantities are 

-2 



where V3 = 2-k denotes the area of the unit three-dimensional sphere. 



s However, the expressions for the regularized action and the conserved quantities will be the same for 
any counterterm choice. 
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Other quantities of interest are the Hawking temperature Tjj = 1/(3 and the area Ah 
of the black hole horizon 



TH = Vhh^ A H = y^r H V 3 . (3.7) 

47T 

The thermodynamics of the EGB black holes can be formulated via the path integral 
approach |3^, ^0|. Following the standard prescription, one computes the classical bulk 
action evaluated on the equations of motion, by replacing the R + \Lqb volume term with 
2(R t t + j(Hj + ^Lcb))- Then one makes use of fl2.14| ) to express the volume integral of 
this quantity as the difference of two boundary integrals. The boundary integral on the 
event horizon is simplified by using the identity (2.13) which provides the following relation 



between the asymptotic parameter W (which fixes the equal-magnitude angular momenta) 
and the horizon data which enter ( |2.15| ): 



H 

A straightforward calculation using the asymptotic expressions (12. 161) shows that the di- 

(E) 

vergencies of the boundary integral at infinity, together with the contributions from D 

and 1^ , are regularized by I c t- As a result, one finds the following finite expression for 
the classical action in terms of the event horizon data and the asymptotic parameters ti 
and W 

t Vz ( ( fTT 2 m , ^ H W \ Vh^ 2 \ 

[^ r * 7^ J + a ^r {hH " irH) ) ■ (3 - 9) 

We have verified that a similar result for I c ; is obtained when using instead the stan- 
dard background subtraction regularization procedure (with a Minkowski spacetime back- 
ground). 

Upon application of the Gibbs-Duhem relation to the partition function, one finds the 
entropy S = (3{E — 2VLhJ) — Id-, which is the sum of one quarter of the event horizon area 
(the Einstein gravity term) plus a GB correction 

S = Se + Sgb, with Se = ^r H y/hH, Sgb = a^y/hii(A- ^). (3.10) 

4G 4G r H 

It is interesting to note that the above expression for the entropy can also be written in 
Wald's form M\ as an integral over the event horizon 

S= IG I d3 ^( 1 + f^)' ( 3 - n ) 

(where h is the determinant of the induced metric on the horizon and R is the event horizon 
curvature) . 
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4. The rotating EGB solutions 



4.1 Non-extremal black holes 
4.1.1 General features 

In the absence of closed form solutions, we relied on numerical methods to solve the EGB 
equations. In this work 9 , we integrated the system of coupled non-linear ordinary differen- 
tial equations with appropriate boundary conditions which follow from ( 2.15| ), ( 2,16| ), by 



using a standard solver J41|. This solver involves a Newton- Raphson method for boundary- 
value ordinary differential equations, equipped with an adaptive mesh selection procedure. 
Typical mesh sizes include 10 2 — 10 3 points. The solutions in this work have a typical 
relative accuracy of 10 -6 . 

In our approach, the input parameters 10 are the GB coupling constant a, the event hori- 
zon radius rjj and the horizon angular velocity £Ih (or equivalently, the equal-magnitude 
angular momenta J through the parameter W). Physical quantities characterizing the 
solutions can then be extracted from the numerical solutions. 

For most of the analysis, we set r# = 1, which does not spoil the generality of the 
results, since it corresponds to fixing the scale of the problem. Also, to simplify the problem, 
we restrict our integration to the region outside the event horizon 11 . 

In constructing rotating EGB black holes, we made use of the existence of the closed 



form solutions ( 2.17 ) and ( |2.20D , and employed them as starting configurations, when 
increasing gradually a or Qh, respectively. Our results clearly show that any such MP 
solution admits generalizations in EGB theory. When starting instead from the EGB static 



black holes (2.20), a corresponding branch of rotating EGB solutions emerges smoothly for 
any value of the event horizon radius r#. 

The profiles of the metric functions of typical EGB black hole solutions are presented 
in Figure la for a static (Qh = 0) and a rotating solution with £Ih = 0.68 (where both 
solutions have the same value of the GB coupling constant). The effects of the GB coupling 
constant are shown in Figure lb, where rotating solutions for two values of a and the same 
horizon angular velocity are shown (here the solution with a = 0.01 is very close to the 
MP configuration). One can see that, apart from a nonzero w(r), the rotation leads to 
non-constant values for h(r)/r 2 , while a nonzero a leads to nontrivial deformations of all 
metric functions. 

These rotating EGB solutions possess also an ergoregion inside of which the observers 
cannot remain stationary, and will move in the direction of the rotation. The ergoregion 
is bounded by the event horizon, located at r = r# and the stationary limit surface, or 
ergosurface, located at r = r e , where the Killing vector d/dt becomes null, i.e. gu = 
—b(r e ) + h(r e )w(r e ) 2 = 0. One can see that the ergosurface does not intersect the horizon. 



9 The numerical methods here are similar to those used in literature to find numerically black hole 
solutions with equal-magnitude angular momenta in Einstein-Maxwell theory J^, EjJ . 

10 We have assumed that Qh > 0, which can always be achieved by t — > — t if necessary. Also, in string 
theory, the GB coefficient a is positive, which is the only case considered here. 

"However, similar to the static case, the GB term is expected to drastically affect the geometry in the 
region inside the horizon of the rotating black holes (see e.g. |J42|). 
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Figure 1: The metric functions f(r), b(r), w(r) and h(r)/r 2 of typical EGB black hole solutions are 
shown for (a) two values of the horizon angular velocity and (b) two values of the Gauss-Bonnet 
coupling constant a. 

The effect of the GB term is to increase the size of the ergoregion. For a fixed event 
horizon radius r h = 1 and = 0.33 one finds e.g. r e (a = 0) ~ 1.059, r e (a = 1) ~ 1.08 
and r e (a = 2) ~ 1.104. 

Another interesting feature concerns the sign of the quantity p e g = f-ff| which, via the 
modified Einstein equations = —jH^ = Tj~, corresponds to a local 'effective energy 
density'. This effective stress tensor, thought of as a kind of matter distribution, in principle 
may violate the weak energy condition. Indeed, this is the case for both black string [13] 
and black ring |I3] solutions in EGB theory. However, we have found that similar to the 
case of the static spherically symmetric black holes ( |2.20| ), the 'effective energy density' is 
a strictly positive quantity for all rotating solutions investigated. This result together with 
those in |[^, 14] leads to the conjecture that in order to have p e g > the solutions should 
possess a spherical topology of the event horizon and be cohomogeneity-1 12 . 

4.1.2 Domain of existence 

Because EGB black holes are found starting both with the MP configurations and with 
the Schwarzschild-GB static solutions, we conclude that, for fixed rjj, rotating EGB black 
holes should exist in a given domain of the (a, 0#) plane. To map out this domain, we 
fixed certain values of a > and increased £Ih gradually. We then obtained numerical 

(c) 

solutions up to a maximal value of the horizon angular velocity , that depends on a. 

The solutions are numerically robust but the integration becomes difficult, when the 
maximal value is approached. However, for any a > 0, no critical phenomenon (like 
e.g. a bifurcation or the approach to a singular point) seems to arise there. 

(c) 

In order to clarify the issue of the limiting solutions for Qh — > £l H for fixed a, the 
study of the event horizon values b'(ru) = b\ and f'{rn) = fi as functions of ru turned 



12 This conjecture could be tested by constructing the d = 5 caged black holes in d — 5 EGB theory (i.e. 
Kaluza-Klein configurations which resemble the Schwarzschild black hole close to the horizon and approach 
asymptotically the four dimensional Minkowski space times a circle). 
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Figure 2: The 'reduced' dimensionless physical quantities (a) horizon area, (b) temperature, (c) 
horizon angular velocity, and (d) ratio between the Gauss-Bonnet and the Einstein gravity contri- 
butions to the total entropy are shown versus the dimensionless squared (equal-magnitude) angular 
momenta for several values of the Gauss-Bonnet coupling constant. The dots represent the data 
points. The spline- interpolated curves have been extrapolated to the extremal endpoints. 



out to be crucial. 

For Qh = 0, we know that b'(rn) = ^ r H/( r H + a )> as can be seen from the explicit 
solution (|2T20| ). Then, increasing £Ih, our numerical results show that b'(rn) and /'(rjy) 
both decrease monotonically. This strongly suggests that they reach the value zero in the 

(c) 

limit £Ih — > & H ■ 

From these observations we conclude, that the families of rotating EGB black holes 

(c) 

terminate at extremal configurations. All relevant quantities remain finite as £lu — > fijy , 
while the Hawking temperature vanishes in the limit. Moreover, we evaluated a number 
of curvature invariants {e.g. the scalar curvature R and the Kretschmann scalar), and our 
extrapolations showed that these stay finite everywhere in that limit. 

This behaviour is thus analogous to that of MP black holes in the extremal limit, where 
= l/y/2rn- We have found this picture for any value of a that we considered. 
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a 

Figure 3: The critical values of Hh, J and E corresponding to extremal limiting solutions are 
shown as a function of the Gauss-Bonnet coupling constant. The dots represent the data points, 
the curves are obtained by spline-interpolation. 

Therefore we conjecture that rotating black holes with equal-magnitude angular momenta 
exist for any value of the GB coupling constant. 

To illustrate these aspects, we show in Figure 2 the behaviour of the 'reduced' area 
of the horizon an, the 'reduced' temperature tu, the 'reduced' horizon angular velocity 
ljh, and the ratio Sgb/Se as functions of the squared 'reduced' (equal-magnitude) angular 
momenta j for several values of the GB coupling constant a. These 'reduced' dimensionless 
quantities are denned as follows 

3 / 3 A H t /27T<7 r- 2 8G ~ „ 2 27vr J 2 fA s 

One can see that the pattern in Einstein gravity 13 is recovered for all values of a considered. 
A nonvanishing GB term decreases, however, the maximal values of the dimensionless 
quantities j 2 , an, and tn, while it decreases the maximal value of loh- 

To determine the domain of existence of the black holes with respect to the GB coupling 
constant a, we exhibit in Figure 3 the values of f2#, J and E corresponding to critical 
solutions as functions of a. (Note, that these quantities are normalized there with respect 
to the corresponding ones for the MP solution). The effect of the GB interaction is clearly 
perceptible on the diagram. One should note that all these quantities have a nontrivial 
dependence on a. 

We conjecture that the domain of existence of rotating EGB black hole is the region 
below the (spline-interpolated) curve for Q H ' (or equivalently below the (spline-interpolated) 
line for the equal-magnitude angular momenta J^). 

13 The following relations hold for the solutions of Einstein gravity 

a H {f) = \{l + y/l=P), tg 0- a ) = 2(l- i^ffHZ ), 4(/).^EZ)„ L (4.2) 
2 3 3 
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Figure 4: The entropy S and mass-energy E are shown vs. the Hawking temperature Th for fixed 
angular momenta J and several values of the Gauss-Bonnet coupling constant (here and in Figures 
5, 6 and 7 we set G = 1). 

To control the quality of these results, we have performed (a number of) these numerical 
calculations in two different gauges, finding excellent agreement for the physical parameters 
of the rotating EGB black holes. Our systematic analysis was limited to a G [0, 3], because 
the GB term is supposed to emerge as a correction to the Einstein lagrangian; however, we 
also found families of solutions exhibiting the same pattern for larger values of a. 

4.1.3 Thermodynamical properties 

Considering the thermodynamics of these solutions, the EGB black holes should satisfy the 
first law of thermodynamics 

dE = T H dS + 2n H dJ. (4.3) 

One may regard the parameters S, J as a complete set of extensive parameters for the mass- 
energy E(S, J) and define the intensive parameters conjugate to them. These quantities 
are the temperature and the angular velocities. Also, for a / 0, these solutions do not 
appear to satisfy any simple Smarr relation 14 . 

In the absence of an exact solution, we attempt here to analyze the thermodynamic 
stability of the rotating EGB solutions based on the available numerical data. Since we 
did not yet explore the full parameter space of solutions, the results below are only partial, 
and one cannot exclude the existence of new features outside the explored domain. 

It is known that different thermodynamic ensembles are not exactly equivalent and 
may not lead to the same conclusions since they correspond to different physical situations 

14 The Smarr relation in Einstein gravity is E = 3(ThS + 2QhJ)/2. The recent work EJJ proposes a 
Smarr relation for static black holes in Lovelock gravity. It would be interesting to extend the formalism 
in [M3| to the case of spinning solutions. 
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(see e.g. |44]], |4q| , [}46f for recent work on the thermodynamics of spinning black objects in 
Einstein gravity). 

In the canonical ensemble, we study black holes holding the temperature Th and the 
angular momenta J fixed, the associated thermodynamic potential being the Helmholz free 
energy F[Th, J] = E — ThS. In this case, the numerical analysis for several values of a 
indicates that the qualitative thermodynamical features of the MP solutions are also shared 
by their EGB generalizations. The configurations near extremality (i.e. with a small enough 
Th) are thermally stable in a canonical ensemble since Cj = Th(OS / 8Th)\j > 0. However, 
there is also a branch of large black holes whose entropy is a decreasing function of Th- 
At the critical point, the specific heat goes through an infinite discontinuity, and a phase 
transition takes place. This is the picture one finds for the MP solution in Einstein gravity, 
in which case the specific heat changes the sign for T H ~ 0.087396/(GJ) 1 / 3 . Interestingly, 
for a given J, the effect of the GB term is to decrease the critical value of the Hawking 
temperature. These features are shown in Figure 4 where we exhibit both the entropy and 
the mass-energy of the solutions as functions of Th for an arbitrary fixed value of J and 
two nonzero values of a as well as for a = 0. 

In the grand canonical ensemble, on the other hand, we keep the temperature and 
the horizon angular velocity fixed. In this case the thermodynamics is obtained from the 
Gibbs potential G[Th, &h] = E — ThS — 2S7# J = I c i/ ft. The first quantity of interest here 
is the specific heat at constant horizon angular velocity Cn = Tn(dS / dTH)\n H > 0. A 
straightforward computation shows that for the Myers-Perry black hole this is a negative 
quantity 

i s/2 + x 2 - 2VT+^(2 + 2x 2 + VTT^) n . + , n H .... 

Co = — — i ; , < 0, with x = — — — . 4.4 

ATTGT H x 2( 1 + 2 .2)3/2( 1 + x 2 + ^/T^2) ' 7rT# 
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Figure 6: The angular momentum J is shown as a function of the horizon angular velocity £Ih 
for a fixed Hawking temperature and several values of the Gauss-Bonnet coupling constant. 



However, not completely unexpected 15 , our numerical results show the existence a branch 
of rotating black holes with Cq > for any value of a > considered. 

To illustrate this feature, we show in Figure 5 the entropy and the mass-energy of 
solutions as functions of the Hawking temperature for several values of the GB coupling 
constant and an arbitrary fixed value of the horizon angular velocity (for better visualiza- 
tion, we used a logarithmic scale for S). We observe that the critical temperature where 
Cr> changes sign decreases with increasing a, while the Einstein gravity picture is recovered 
for large black holes, where Cq, < 0. 

Another response function of interest is the isothermal permittivity €t h = (dJ/dQii)\T H - 
The Einstein gravity result is 

V 3 5 + 6x 2 - \/l + x 2 (5 + 3x 2 ) , n H -x 

er„ = -r— — - A ; ; — , where x = — — — , 4.5) 

lH W^GT^ x 2 (l + x 2 ) 3 / 2 (l + VTT^ 2 ) 2 kTh K ' 



and thus ex H > for Qh/ttTh < 0.77272. Our results for a = 0.1 and a = 1 show that 
changes sign also for EGB solutions and that, for a given Th, the solutions with a small 
{In have a positive isothermal permittivity (see Figure 6). Thermodynamic stability in the 
grand-canonical ensemble requires that the specific heat at constant angular momentum 
Cj, the isothermal permittivity ct h as well as the specific heat at constant horizon angular 
velocity are positive. Although a systematic study of the domain of thermodynamic 
stability of the EGB solutions is beyond the purposes of this paper, our results show the 
existence of solutions which fulfill these conditions. The stability conditions hold for a 
set of rotating solutions emerging from small static black holes (that have positive specific 
heat). 

We expect that this result is relevant for the properties of the corresponding d = 6 
rotating black string solutions in EGB theory (note that such solutions cannot be found 



3 We recall the existence of a branch of small static EGB black holes that is thermodynamically stable 
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by simply uplifting the black holes in this paper). For example, the results in [47] indicate 
that the Gregory-Laflamme instability Ej| persists up to extremality for all d = 6 Einstein 
gravity black strings with equal magnitude angular momenta. Based on the Gubser-Mitra 
conjecture [49] that correlates the dynamical and the thermodynamical stability, we con- 
jecture that in a certain region of parameter space there are spinning EGB solutions which 
do not possess a Gregory-Laflamme instability. Also, not completely unrelated, it would be 
interesting to study the thermodynamic stability of these solutions from the perturbative 



corrections to the gravitational partition function (see 1 5C ] for a recent similar study of the 
Kerr-AdS black holes). 

We conclude that the presence of a GB term in the lagrangian affects the thermody- 
namical properties of the solutions and allows for thermodynamically stable solutions, both 
in the macrocanonical and the canonical ensemble. 

4.2 Extremal solutions: near horizon geometry and the entropy function 

Returning to the issue of the extremal solutions, we recall that the numerical integration in 
the neighbourhood of extremal black holes is very difficult, as noted also in other cases. The 
near horizon expansion of extremal configurations is still given by ([TL5|), with fx = bi = 0, 



and thus they satisfy a different set of boundary conditions at r = rjj than the nonextremal 
configurations obtained in this work. Therefore finding such solutions explicitly is beyond 
the scope of this paper 16 . Most likely, this should require a different parametrization for 
the metric, better suited to the extremal case. 

However, we argue that the existence of the EGB generalizations of the extremal MP 
black holes for any a > is strongly suggested by the existence of an exact solution 
describing a rotating squashed AdS2 x S s spacetime. This solution would describe the 
neighbourhood of the event horizon of an extremal black hole. (The far field expression 
of the extremal solution is still given by ( |2.16| ), with a single essential parameter in the 
expansion there.) 

Therefore we consider the following metric form (see the generic ansatz ( |2.9[ )) in coro- 
tating coordinates 



ds 2 = v 1 (^-- r 2 dt 2 ) + ^-{af + a 2 ) + ^(a 3 + 2krdt) 2 (4.6) 



{i.e. for b(r) = v\r 2 , f(r) = r 2 /vi, g(r) = v^, h(r) = V2V3, w(r) = kr within the 
parametrization in this paper), such that the horizon is located 17 at r = 0. This ge- 
ometry describes a fibration of AdS% over the homogeneously squashed 5 3 with symmetry 



group 50(2,1) x 5*7(2) x U{1) |3 



The parameters Vi , k satisfy a set of algebraic relations which result from the EGB 
equations. In what follows we choose to determine them by using the formalism proposed 
in 123] , thus by extremizing an entropy function. This allows us also to compute the entropy 
of these black holes and to show that the solutions exhibit attractor behaviour. 

16 Note, however, that we could construct with relatively good accuracy near- extremal black holes. 
17 This position of the horizon can always be obtained by taking r — > r — ru- 
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Therefore let us denote by f(k, v) the lagrangian density y/—gC evaluated for the near 
horizon geometry fl4.6|) and integrated over the angular coordinates, 

J &d4tyy/=t{R+jL GB ). (4.7) 

The metric field equations in the near horizon geometry ( [4.6| ) now correspond to = 
J, ^j- = 0, with J the angular momenta of the solutions. 



f(k,v) = / dOd^y/^C 



Then, following [24|, we define the entropy function by taking the Legendre transform 
of the above integral with respect to the parameter k, 



£(J,k,v) = 2ir(Jk- f(k,v)). 



(4.8) 



It follows as a consequence of the equations of motion that the constants k,v are solutions 
of the equations 



^ = 0, ^ = 0. 

dk ' dv; 



(4.9) 



Then, the entropy associated with the black hole is given by S extre mai = £{J, k, v) evaluated 
at the extremum ( [4.9[) . Further details on this formalism and explicit examples are given 



e.g. in H, |2|, @. 

For the metric ansatz (|4.6|), a straightforward calculation gives 



Jk _ vr^p^ 2 ^ _ ^ 2 _ _ A ^ _ a{k 2 V2V ^ _ 4) _ 4ui(v3 _ 4) ) (410) 
IbGvi 



£(J,k,v) = 2vr 

such that the explicit form of the equations fl4.9|) is 
dS 



dv\ 
88 
dv 2 
d£ 
dv 3 



-I6v 2 + 4vfv 3 + k 2 v\v 3 + ak 2 v 2 v 3 {4 - 3v 3 ) = 0, (4.11) 
=>- -16vf + 12v\v 2 + Av\v 3 - 5k 2 v 2 v 3 + a(-4vi(v 3 - 4) + 3k 2 v 2 v 3 (3v 3 - 4)) = 0, 
=>- -16v 2 + Aviv 2 + 12w^ 3 - 3A; 2 w^ 3 + a(-4vi{3v 3 - 4) + 3fc 2 w 2 i'3(5w3 - 4)) = 0, 



and 



d£ n T irk(v 2 v 3 ) 3 / 2 

^ = 0^J=^-^(, 2 + a(4-3, 3 )). 

In Einstein gravity the solution has a simple form in terms of J 

2/3 



(4.12) 



1 fGJ\ 2/3 ( V2Gf 



, v 3 = 2, k=-, 



(4.13) 



and thus the known result £ = S, 



extremal 



7T J is recovered. 
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In the limit of small a, one can treat the GB term as a perturbation. In second order 



in a, one finds the following solution of the system (4.11), (4.12) in terms of the global 
charge J 



1/GJ\ 2/3 a ( tt \ 3 / 2 a 2 



fV2GjV /3 . 4 . /2;rY i2 -_>. 



v 3 = 2 - 2 



while 



GJJ 



a 



/ vr \4/3 2 1 / 3 16a 2 
\Gl) 3 ' k 



1 / x \ 2 ' 3 



a / 7T \ 4 /3 a 2 



2 + V^2GjJ 2 + Uj; 2 2 / 3 12 ' 



£ = vrJ + 3 



J^\ 1/3 

2G 2 y 



a: 



V 4G 4 J 



" ^a 2 
2 



(4.15) 



Therefore the approximate solutions exhibit a complicated behaviour in terms of a and J 
and a perturbative approach may be misleading. 

Unfortunately, for a nonzero GB term in the action, it seems that the only possibility is 
to express the nonperturbative solution of the system ( 4.1 1| ) , ( |4.12j ) in terms of the relative 
squashing parameter v 3 , with 



J - 



and 



(v2 + (4 — 3^3)0) (3v2 — (i>3 — 4)a) 
2(v 3 - 4)(3v 2 + (8 - 6v 3 )a) 

^^^(4-v 3 )(v 2 + a(A-3v 3 )), 



8GJvi 



tt(v 2 v 3 ) 3 / 2 (v 2 + (4 - 3v 3 )a) ' 
while the radius of the round S 2 sphere in the line element (|4.6| ) is given by 18 



(4.16) 
(4.17) 

(4.18) 



a 



>'2 



v 3 -2 



2v 3 - 7v 3 + 4 - J5v% - 34^ + 73w| - 56v 3 + 16 



(4.19) 



Also, we notice that the relations ( [4.16j )- (l4.19| ) are invariant under the scaling 

v\ ->■ \v\, v 2 ->■ Xv 2 , £ — > A 3//2 £, J — > A 3 / 2 J, k —f k, and a — > Aa, 

which shows that the solutions exist for any a > 0. 

Inserting these expressions into Eq. ( 4.10|) we obtain for the entropy function of the 
extremal black hole: 



7T 



£ = S extr emal = 7^\/v 2 V 3 (v 2 - (v 3 - 4)«) , 



(4.20) 



(with ^2(^3) as implied by ( 4.19j )) 



'Here we restrict to the physical solution which recovers the general relativity limit asa^O. 
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Figure 7: The dimensionless quantities £/a 3 / 2 , k, V\/a, v-ija and V3 are shown as functions of 
J/a 3 / 2 (for a better visualisation, we employ here logarithmic scales). 



As a check, we note that the result (4.20) agrees with Wald's form ( |3.11 ) evaluated 
for the near horizon geometry (|4.6|). We also note that, in principle, £ can be expressed in 
terms of the conserved charge J by inverting relation ( [4.17 ), 

We conclude that a nonzero a may substantially affect the near horizon geometry of an 
extremal black hole. For example, the allowed range 19 for the relative squashing parameter 
V3 is < ^3 < 2 (i.e. the GB term reduces the relative squashing of the solutions), with 



Vl 



- + 

A T • • • 



v 2 



cur, 
2 



+ 



J 



2G + ■■■ 



k 



± L f 



V27T 2 

G 



(^3) 3/2 + • • • , 



as — > 0, and 



2-v 3 



4a 
2-v 3 



+ .., J 



4V2n 
G 



\3/2 
2=^) + " 



k 



i+ £ 



4V2n 2 
G 



n 3/2 

2-V3 J 



as V3 — > 2. In Figure 7 we exhibit a number of relevant dimensionless quantities as functions 
of the scaled angular momenta J/a 3//2 (where for better visualisation, we use logarithmic 
scales). We observe that the ratio £/J is no longer constant in EGB gravity. 

However, finding local solutions in the vicinity of the horizon does not guarantee the 
existence of global asymptotically flat solutions. Further progress in this direction seems 
to require an explicit construction of the bulk extremal black hole solutions. For example, 
this would also allow to construct the E(J) diagram for such configurations. 

We close this Section by remarking that the study of these AdS2 x S 3 solutions in 
EGB theory is interesting in yet another context. In ref. f)5|| it has been proposed that 
the near horizon geometry of an extremal Kerr black hole is holographically dual to a 2- 
dimensional chiral conformal field theory (CFT). This correspondence has been extended 
to various other examples of extremal spinning black holes in d > 4 dimensions, including 



9 This results from the physical condition that Vi, J 2 are strictly positive quantities. 



- 20 - 



configurations with matter fields. These studies are based on the universality character 
of the near horizon geometry of extremal black holes. It would be interesting to consider 
also the case of such solutions with GB corrections, the configurations in this work being 
perhaps the simplest relevant example. In a Kerr/CFT context, the entropy formula ( 4.2C] ) 



should be recovered by computing the central charge of a certain two-dimensional conformal 
algebra 20 . 

We hope to return to the study of extremal black holes in EGB theory in future work. 
5. Further remarks 

The main purpose of this paper was to present evidence for the existence of rotating black 
holes in d = 4 + 1 EGB theory. Representing generalizations of a particular class of MP 
black holes, the considered configurations possess a regular horizon of spherical topology 
and two equal-magnitude angular momenta. Our results indicate that the inclusion of a 
GB term in the action does not affect most of the qualitative features of the solutions. 
However, the presence of a GB term in the action has a tendency to stabilize the rotating 
black holes, leading to a branch of solutions with a positive specific heat at constant angular 
velocity at the horizon which does not exist for the MP solutions. 

Also, analogous to the case of Einstein gravity, when the horizon angular velocity is 
increased, the black hole solutions reach a limiting extremal black hole with a regular 
horizon. Although we did not attempt to construct these extremal black holes, we gave 
further support for their existence in Section 4.2 by finding an exact AdS2 x S* 3 rotating 
solution in EGB theory. This solution would describe the neighbourhood of the event 
horizon of an extremal EGB black hole. 

A natural question that arose during our study was how to determine the various con- 
served quantities and the total action of the solutions. One particularly powerful approach 
to this problem is given by the counterterms method. In Section 3 we showed how to gener- 



alize the Einstein gravity counterterms in 1 2£ , [3(J, 31 1 by including the effects of a GB term 



in the bulk action. Although the counterterm method gives results that are equivalent to 
those obtained using the background subtraction method, we employed it here because it 
appears to be a more general technique than background subtraction. Moreover, it is inter- 
esting to explore the range of problems to which it applies, in particular for configurations 
with higher order curvature corrections in the gravity action. 

The solutions obtained in this paper may provide a fertile ground for the further study 
of rotating configurations in EGB theory. For example, their generalization to include 
the effects of an electromagnetic field is straightforward. Also, in principle, by using the 
same techniques, there should be no difficulty to construct similar solutions in d = 2N + 1 
dimensions with N > 2 equal-magnitude angular momenta. 

Another interesting direction to consider in future work consists in finding the EGB 
generalizations of the d = 5 MP black holes with nonequal angular momenta, in particular 



Note that one can associate a temperature T = l/2nk to the near horizon geometry (4.6). As can be 
seen from Figure 7, in EGB theory this temperature is no longer constant T = presenting a a nontrival 
dependence on the dimensionless ratio GJ/q 3//2 . 
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the case with rotation in a single plane. Concerning the latter case, it would be interesting 
to see how the GB term would affect the properties of the solutions close to the extremal 
limit. Different from the solution ( |2, 17 ), in this limit the Einstein gravity solution corre- 



sponds to a naked singularity. One might speculate that the higher derivative terms in the 
action might smoothen this singularity and lead to a physically reasonable solution. These 
aspects together with the investigation of the effects of the GB term on balanced black 
rings and thus the phase diagram are presently under investigation. 
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